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I. INTRODUCTION 

In a recent series of articles , we explored a new approach to the gravitational 

effects of vacuum energy density. This approach starts from a conserved microscopic variable 
g, whose statics and dynamics are studied on macroscopic scales. 

The precise nature of q is uncertain for the moment, but we have presented at least 
one concrete example in terms of a four-form field F. This F field could be a part of the 
(unknown) fundamental theory of elementary particle physics with an energy scale given by 

^Planck -10^9 GeV. 

In this article, we do not contemplate ultrahigh energies but propose an explicit realization 
of q from well-established physics, namely, quantum chromodynamics (QCD) with an energy 
scale of the order of 1 GeV. That is, we find that q can be identified as a particular gluon 
condensate in the nonperturbative QCD vacuum. Neglecting QCD effects, an F-type field [2I 
may still be required to reduce the macroscopic vacuum energy density from a natural value 
of the order of (-Epianck)^ to a value which is essentially zero . 

With our general understanding of g-theory, we can then investigate the gravitational 
effects of the QCD vacuum. Most importantly, we find that the dynamics of this gluon con- 
densate in the nonequilibrium context of the expanding Universe may result in a nonzero 
limiting value of the vacuum energy density. This remnant vacuum energy density may 
correspond to the inferred cosmological constant responsible for the observed "cosmic accel- 
eration" (cf. Ref. {4] and other references therein). 

In order to be clear about the terminology, we consider a time-dependent gravitating vac- 
uum energy density Pvac(^) in an expanding Friedmann-Robert son- Walker (FRW) universe 
[typically, Pvac(^) decreases with cosmic time t] and define the cosmological constant A as 
the remnant vacuum energy density in the limit of large cosmic times, A = limt_^oo Pvac(i)- 
This also implies that, while the vacuum energy density pvac may have changed with time, 
the equation-of-state parameter w^sjc = -Pvac/Pvac has kept the value —1, at least for the type 
of theories considered here. The possibility of having time-dependent Pvac(^) and constant 
"i^vac = — 1 may be an important lesson for observational cosmology. 

Throughout this article, natural units are used with c = h = 1 and Newton's gravitational 
constant Gn is shown explicitly. 
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II. GLUON CONDENSATE 



The underlying theory of the strong interactions is nowadays beheved to be given by a 
particular non-Abelian gauge field theory called quantum chromodynamics; see, e.g., Ref. [sl 
and other references therein. The non-Abelian gauge group is SU (Ac) and the perturbative 
particle content of QCD is given by Af fiavors of quarks and A^ — 1 = 8 types of gluons, for 
Ac = 3 colors of each quark fiavor. The nonperturbative particle content of QCD is given by 
the genuine asymptotic states, that is, the baryons, the mesons, and possibly the glueballs. 

The crucial object for our discussion is the Yang-Mills field strength, defined as 

GV(^) = d,A\{x) - d.A\{x) + r'^A''^{x)A%{x) , (2.1) 

with spacetime indices fi, v ranging over to 3, Lie-algebra indices a, 6, c taking values from 
1 to N"^ — 1, repeated Lie- algebra indices 6, c being summed over, and structure constants 
jafec corresponding to the Lie algebra su(Ac). Note that the gauge coupling constant g has 
been absorbed in the gauge potential, so that A"- ^ = 0{g^) for an instanton configuration 
and A"-^ = 0{g^) for a perturbative configuration with a few gluons. 

Consider, now, the QCD gluon condensate (see Ref. 6|| and other references therein). It 
follows directly from gauge invariance that 

(0|GV(x)|0) = 0, (2.2) 

where the vacuum expectation value can, for example, be obtained from a Euclidean path 
integral. Equation (12.21) for the vacuum expectation value of a Yang-Mills field of local 
support is a special case of Elitzur's theorem j^], which relies on the gauge noninvariance of 
the Yang-Mills field strength and the gauge invariance of the path integral and the vacuum 
state. As a further clarification of (12. 2p . we state our explicit assumption that the so-called 
Savvidy vacuum |8i] is not realized, i.e., that the vacuum expectation value of the average 
color magnetic field is zero. 

The vacuum expectation value of the quadratic expression can, however, be nonzero: 

(0| ^ ^V(^) |0) = ^(^) (^?MP(^) 9uAx) - 9,aix) g,,{x)) , (2.3) 

again with an implicit sum over the repeated Lie-algebra index a. The explicit realization 
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of the vacuum variable q from Ref. l| is then 



g(x) = (0|^G"^^^(x)GV(^)|0), (2.4) 



which, with the chosen numerical factor, is precisely equal to the Shifman-Vainshtein- 
Zakharov condensate as determined from charmonium data [q ~ 10^^ GeV^ ) . On the 
theoretical side, note that the vacuum expectation value (12. 4p is a properly renormalized 
quantity, which can, for example, be obtained from a Euclidean path integral calculated 
in the dilute-instanton-gas approximation (see Sec. 6.7 of Ref. 0, (a)] and other references 
therein) . 

The experimental value for q is positive, even though expression (12.41) is not positive 
definite for a Lorentzian spacetime metric. However, q is manifestly positive definite for a 
Euclidean spacetime metric, which is anyway needed to make sense of the path integrals 
for instanton-type calculations. Henceforth, we consider the vacuum variable q to be non- 
negative. 

Before we start our discussion of the gravitational effects of the gluon condensate, we can 
already mention a side-product of our investigation, namely, that the gluon condensate has 
a new characteristic, the compressibility x- This will be mentioned briefly in Sees. |V]and|VT] 
and Appendix lAl while the general discussion of vacuum compressibility has been presented 
in Ref. 



III. COSMOLOGICAL TERM FOR GRAVITY 



The goal of the present article is to explore certain gravitational effects of the QCD gluon 
condensate over spacetime volumes very much larger than those corresponding to the typical 
scales of QCD, £qcd = ctqcd ~ Ifm = IQ-^^m ^ he/ {2m MeV). In the spirit of Ref. Q, 
we consider the following coarse-grained effective action: 

Ses[g, q] = 5'grav[^] + 5vac[^, ^] = j ^"^^ ( ^q^q^ ^Id] + e(^/)) , (3.1) 

where the pure-gravity action Sgrav is given by the standard Einstein-Hilbert term with the 
Ricci curvature scalar R = R[g] and the vacuum-energy- density action S'vac is determined 
by a general function e{q) of the vacuum variable q. Here, q is realized as the vacuum 
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expectation value (12.40 but now averaged over a (Euclidean) spacetime volume of the order 
of (Ifm)^ 

The energy-momentum tensor obtained by variation over g'^'^ is then given by 

2 SS^^c . . ^ de{q) 6q 

= -7=9 = '^'^ - ' ^ 5^ ■ 

Using (12. 3p and (12. 4p . one obtains 

^ G'^^.G^ \^?'"^ = ^g<7M., (3.3) 



5^/.^ \47r2 w 1 ^ 2 
where the brackets of the expression in the middle denote both the vacuum expectation 
value and an average over a (Euclidean) spacetime volume of the order of (Ifm)"^ [the 
same expression can also be written in terms of the effective fields (14. ip introduced below]. 
As a result, (13. 2p produces a cosmo logical-constant-type energy- momentum tensor for the 
Einstein field equation, 

^t^Aq) = pY!^{.q) Qfiu , (3.4a) 

Pvac(g) = e(g)-g^, (3.4b) 

where the expression for Pvac{q) has precisely the structure argued on thermodynamical 
grounds in Ref. The term (I3.4ap in a cosmological context corresponds to a cosmic fluid 
with equation-of-state parameter w^ac = — 1- 

At this moment, it may be instructive to comment on the difference between our approach 
and the one of nonlinear electrodynamics as discussed in, e.g., Refs. [9|, |10|]. In our approach, 
the average energies of the magnetic and electric field fluctuations in the vacuum are related 
by (|Bp) = — (|Ep) = g/4, as follows from (12. 3p . (The negative value of the average 
energy of the electric field is obtained after renormalization of the divergent energy of the 
quantum fluctuations.) But, in the approach of Refs. [o], [lo|, both quantities are non- 
negative, ( |Bp ) > and ( |Ep ) > 0, so that generically the energy-momentum tensor from 
the electromagnetic field does not correspond to a cosmological-constant-type term (I3.4ap . 



IV. EQUATION FOR q 



The equation of motion for q can be obtained by averaging the effective Yang-Mills 
equation. We proceed by the introduction of a "master gauge field" (denoted by a bar), 
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with the following properties 



G\„ = d,A\-d.A\ + r'^A\A\ 



(4.1a) 
(4.1b) 



Physically, the idea is that the classical master field describes the gluon condensate and 
allows for variations over spatial and temporal scales which are large compared to the mi- 
croscopic scales of QCD. Theoretically, such a master field is known to exist in the large-iVc 
limit; cf. Ref. ll|. In a follow-up article, the gradient-expansion method will be used, which 



allows us to get more explicit results [in this method, the vacuum order parameter q{x) is 
considered to be a slow (hydrodynamic) variable with a length scale of inhomogeneities large 
compared to the QCD length scale]. 

Now, start from the variation of (13.11) with respect to A^{x) = A^{x) , where the T" 
are the anti-Hermitian generators of the Lie algebra su{Nc). The variational principle then 
gives the following field equation: 



(4.2) 



where q appearing in the function e'(g) = de/dq stands for the expression (I4.1al) and 
denotes the covariant derivative with respect to general coordinate transformations [using 
the standard affine connection r^^,(a;)] and non-Abelian gauge transformations [using the 
master gauge field A"'^{x)]. 

Next, contract (14. 2 p with G^i/ = G^'^^^T"-, multiply by l/(47r^), and take the trace (with 
normalization factor —2): 



-2tr 







(4.3) 



Using 



2tr [G..Gn = G\,G''''' = q^ 5/ 



(4.4) 



one obtains 



— / de{q) \ ^ ^ dejq) 
\ dq J dq 



2tr 



Gku D^G 



47]-2 



(4.5) 
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At this moment, we can proceed in two directions. The first direction assumes that the 
physical situation is such that the master field satisfies the standard Yang-Mills equation, 
D^j^G^" = 0. Then, the following result holds: 

ti [G^^D^G'^''] = 0, (4.6) 

which nullifies the right-hand side of (14. 5p . 

Equation (14.61) for the case of Minkowski spacetime can also be argued as follows. Take for 
granted that the nonperturbative QCD vacuum over Minkowski spacetime does not break 
spacetime translation invariance and also does not break any of the discrete symmetries of 
charge conjugation (C), parity refiection (P), or time- reversal (T). (The implicit assumption 
is that the so-called 6 parameter vanishes; cf. Ref. 5|].) Then, it follows that the k = com- 
ponent and the k, = 1,2,3 components of the left-hand side of (14.61) vanish by, respectively, 
the T and P invariance of the Minkowski-spacetime QCD vacuum. 

As g is a gauge-invariant scalar, the covariant derivative on the left-hand side of (14. 5 p 
equals the standard gradient and the solution of (14. 5 p using (14. 6 p is simply 

(4.7) 

aq 

where fi is an integration constant. [It will be shown in a forthcoming publication that 
(14. 7p also follows from the gradient expansion up to the first-order (linear) term in ^^g.] 
Result (14. 7p demonstrates that the density q of the gluon condensate is a conserved quantity 
and that fi from (14.70 plays the role of the corresponding chemical potential. The physical 
situation corresponds, therefore, to that of an equilibrium state of the vacuum. 

The second direction considers a physical situation with additional higher-derivative 
terms contributing to the equation of motion for the master field, so that D^G^'^ need 
not vanish in general. The equation of motion (14.51) can then be rewritten as 



-2tr [l/{4n') G.^D.G"'] 
d^-jl = -4/1 ^7 _ / , (4.8a) 



2tr[l/(47r2)G,,,Gn 

where 7i(x) is an effective gauge-invariant scalar field and the denominator on the right-hand 
side of (I4.8ap is precisely equal to q according to (14. 4p with its spacetime indices contracted. 



The dynamical equation fl4.8ap will be used in Sec. |VT] and Appendix |X] to estimate the 
remnant vacuum energy density for the present (T-noninvariant) Universe. 



V. EFFECTIVE POTENTIAL FOR q 



In the simplest approach, the effective potential for q is determined by asymptotic free- 
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12! and the conformal anomaly [13| 



and the conformal anomaly [13| evaluated at one loop: 
q 

= eo + &ig In — , (5.1a) 

<lc 

with number of colors N^. = ^ and number of light-quark flavors N{ = 2 for QCD at low 
energies [recall that q as defined by ( 12. 4p contains an explicit factor l/(47r^)]. 

With this choice for the effective potential e(g), (14.71) gives the following expressions 
for the gluon-condensate charge q, the macroscopic vacuum energy density Pvac, and the 
energy-momentum-tensor trace T,^ as a function of the chemical potential /i: 

q{n) = qc exp [fi/bi - l) , (5.2a) 

Pvac(/i) = e(g(/u)) - l^qil^) =eo-bi g(/i) , (5.2b) 

T/(/i) = 4eo-46ig(/i). (5.2c) 

The second term on the rig ht-hand side of fl5.2cp corresponds to the conformal anomaly, as 



discussed in, e.g., Ref. [ij], where /x is the chemical potential of baryons and reflects the 
conservation of baryonic charge. 

Here, fi is the chemical potential that characterizes the vacuum state and reflects con- 
servation of the vacuum charge q. Moreover, fi becomes a "running coupling constant," 

/i = 61 (1 + ln(g/g,)) , (5.3) 

as follows from (I5.2al) . The gluonic vacuum is stable, since the vacuum compressibility |l| 
is positive for 61 > and g > 0: 



x = [q^ 



dq^ 

as will be discussed further in the next section 
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{hq)-\ (5.4) 



VI. A FROM A SELF-SUSTAINED GLUONIC VACUUM 



Given that q from fl2.4p and hi from fIS.lbl) are non-negative for low-energy QCD, the 
vacuum energy density Pvac(/u) in (15.2bp can be nulhfied if cq > 0. In this case, the self- 



sustained vacuum is given by 

Po = ^1 (^1 - In&i + In (eo/gc)) , (6.1a) 
g(/io) = qo = eo/&i , (6.1b) 

Pvac(/io) = --Pvac(po) = 0, (6.1c) 

r/(/xo) = 0, (6.1d) 
where the result for the vacuum pressure in fl6.1cp follows from the general energy-momentum 



tensor fl3.4al) . Recall that a self-sustained vacuum jl| can exist as an equilibrium state 
at zero external pressure Pext, with pressure equilibrium giving P^ac = -Pcxt = 0. The 
particular gluon-condensate vacuum discussed here has a vacuum compressibility (15. 4p given 
by xo = x(Q'o) = 1/^0 > 0, according to fl6.1bp . 



For the case eo < and with g > 0, the energy density can only be nullified if 

hi < 0, which holds for an Abelian gauge field theory such as QED. The vacuum would, 
however, be unstable, since the vacuum compressibility (15.41) would be negative for negative 
hi. In addition, it is far from obvious that a nonzero vacuum expectation value (12.40 for q 
can arise in an Abelian gauge field theory. In short, a stable self-sustained vacuum can be 
reahzed by a non- Abelian gauge field theory with eo > but not by an Abelian gauge field 
theory. ^ 

The quantities g(/io) and qc are determined by the characteristic QCD energy scale Aqcd 
from the asymptotic-freedom behavior [l2| of the S'[/(3) gauge coupling constant, 

g(/io) ~ gc ~ A'qcd - (200 MeV) \ (6.2) 
with eo ~ hi Aq^j-, from fl6.1bp . Still, the macroscopic energy density of the self-sustained 



^ The quantity eo may, in principle, come as the response of (or reaction from) the deep vacuum at the 
Planck-energy scale, which is slightly adjusted to compensate the energy of the gluon condensate, as 
discussed in Sec. II of Ref. [l,] for the case of a scalar condensate. For the present discussion, cq is simply 
assumed to be positive. As mentioned above, cq then corresponds to the inverse vacuum compressibility. 
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vacuum that enters the Einstein equation as a cosmological constant is not given by e(go) ~ 
eo = 0(10^^ eV^) but is strictly zero, A = Pvac(/io) = 0, according to (1^ . (HTj) . and (KTch . 

A nonzero value of A = p^adf^) may appear for a perturbed vacuum with fi 7^ fiQ. 
Specifically, the vacuum energy density induced by the expansion of the Universe can be 
expected to be nonzero (cf. Sec. [IV|l . Based on the heuristic discussion in Appendix |X1 we 
suggest the following behavior: 



p.a,~/|i/|A^QCD, (6.3) 

with Hubble parameter H = {da/dt)/a > for an expanding universe and a numerical 
factor / > 0. It has indeed been argued 15j on general grounds that the linear H term of 
fl6.3p may arise from the nonperturbative QCD interactions that anomalously break the scale 
invariance of the massless classical theory. The potential importance of the QCD vacuum 
for cosmological horizons has also been emphasized in Ref. 

According to our present understanding (see, e.g., Ref. the Universe evolved from 
an early radiation/matter-dominated phase [H{t) ~ 1/t] to a late vacuum-dominated phase 
[H(t) ~ const.]. The crossover will be discussed further in the next section, but, here, only 
the asymptotic behavior {t —>■ 00) will be considered. 

For a stationary de-Sitter universe, result fl6.3p can be written as 

A = /ifdesA^CD, (6.4) 

with Hdes > the Hubble constant (time-independent Hubble parameter) of de-Sitter space- 
time and neglecting higher-order terms such as -f^des ^qcd- addition, the standard Fried- 
mann equation gives for a de-Sitter universe 

^^des = (8vr/3) A/E^^_, , (6.5) 



with -Epianck = v^^c^/Gn ~ 1-22 X 10^^ eV. Eliminating H^es from the last two equations, 
one obtains the following estimate of the cosmological constant (remnant vacuum energy 
density) : 

A=(87r/3)f A6QCD/i?2^,„,k, (6-6) 
where the numerical constant remains to be determined. 
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As the QCD scale parameter Aqcd of the SU{3) gauge couphng constant is 
renormahzation-scheme dependent, it may be more appropriate, conceptually, to give the 
cosmological constant in terms of a directly measurable quantity. Specifically, we take the 
string tension Kqqd = 1/ (27ra') ^ (400 Me V)^ from the measured Reg ge s lope a' of meson 
resonances or from numerical calculations of lattice gauge theory 17 1 combined with 
other experimental data to fix the absolute length scale. Setting Aq^jj) ^ Kqcb/4: and 

^ {24/ it) k\ in (16. 6p . the final expression for the cosmological constant reads 
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A ^ A-ac„/^^?,.„. « (3 X 10-ev)' (^) (^^^ J . (6.7) 

where the numerical constant k\ remains to be determined (the experimental results to be 
discussed shortly suggest a value of the order of 10~^). Result (16. 7p can also be written as 
A ~ (Gn/c^) -^qcd/^; order to emphasize that the result relies only on classical general 
relativity and quantum chromodynamics (the string tension Kqqd has the dimension of 
energy over length). 

The suggestion, then, is that the vacuum of the presently observed Universe is not relaxing 
to the absolute equilibrium state (16.11) but to the de-Sitter equilibrium state with nonzero 
cosmological constant (16.61) or equivalently (16.71) . 

Since the proton mass nip ^ 938 MeV is now known to come mostly from the gluon 
dynamics, nip ~ Aqcd, estimate (16. 6p corresponds to Zeldovich's original suggestion [18] for 
the cosmological constant in terms of the proton mass, Azeidovich ~ "^p/-^pianck- Numerically, 
one has ^ ^qcd Zeldovich's expression gives a value for A several orders of magni- 
tude larger than (16. 6p . which is closer to the observed value but still somewhat too large for 
/ = !• 

The numerical a gree ment between the theoretical estimate (16.60 or (16. 7p and the ex- 



perimental value 



20 



21 1 of approximately (2meV)^ is improved by having a reduction 



factor k\ = 0(lO~^) in (16. 7p . The corresponding factor / = 0(lO^'^) in (16. 6p traces back 
to (16. 3p and depends on the evolution of the gluon condensate as the Universe cools from 
T ^ 200 MeV to the present temperature T ^ 3K; see Appendix lA] for further details. 
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VII. OTHER COSMOLOGICAL CONSTANT PROBLEMS 



In the previous section, we have made an 
cosmological constant problems (cf. Refs. [3, 



attempt to use QCD for the following three 
nd other references therein): 



(i) why is the cosmological constant A not catastrophically large? 

(ii) why does not A vanish exactly? 

(iii) what physical mechanism sets the scale of A? 

From the g-theory approach to QCD, we have found that the Universe asymptotically 
approaches a stationary de-Sitter phase with a cosmological constant A given by fl6.6p or 
equivalently (16 ■71) . which suggests a partial solution to the above three problems.^ 

However, not all cosmological constant problems have been addressed, let alone solved 
completely. There remain, for example, the following two questions: 

(iv) at which moment in time, t = tcross, does the vacuum energy density start to dominate 
over the cold-dark-matter energy density? 

(v) why do galaxies and stars exist at times relatively close to tcross? 

Question (iv) perhaps has a simple answer in our approach. The crossover from the 
cold-dark-matter-dominated Universe to the gluon-condensate-dominated Universe occurs 
when the cold-dark-matter energy density drops below the vacuum energy density. For a 
fiat matter-dominated FRW universe with Hubble expansion parameter H ^ (2/3) the 
cold-dark-matter energy density evolves as Pcdm(^) ~ (S/Svr) (4/9) Ep^^^^^^^f t"^ , with numerical 
factors of order unity displayed.'^ Asymptotically {t ^ tcmss), the constant vacuum energy 
density is given by (16. 6p . The resulting crossover time can, therefore, be estimated as 



Across ~ (47r) / ^ -Epianck/^QCD • 



(7.1) 



^ As mentioned in Sec.IH non-QCD contributions to the vacuum energy density are perhaps canceled by the 
self-adjustment of another q-type field such as the 4-form field F considered in Ref. [2] or by an entirely 
different mechanism. 

^ In first approximation, the energy transfer from vacuum to cold dark matter can be neglected for t <C Across i 
as long as the vacuum energy density is given by (|6.3p also for a time-dependent Hubble parameter H. 
The question remains as to the precise nature of the energy exchange between vacuum and matter [19 1. 



12 



In terms of the string tension Kqcb ~ 4 Aqqq and the numerical constant ~ (7r/24) 
the crossover time (17.11) becomes 

«o™ « (6. ^.r'' ^k.. A-^g 2 X 10" = [^) " ( " . (7.2) 

where a value for /ca of the order of 10~^ is indicated by the comparison of (16.71) with the 
measured vacuum energy density, as discussed in Sec. IVll 

The value for tcross from (17. 2p is of the order of and even just under the observed value 
for the age of the present Universe, to ~ 14Gyr ^ 4 x 10^'' s, as determined from the data 



compiled in Refs. 
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2l|. The corresponding redshift -Across = 0(1) agrees with the results 
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indicated by deep supernovae observations, such as those reported in Ref. 22l].^ 

Question (v) remains unanswered for the moment, but the answer could also be related 
to QCD, possibly via the mass and baryon number of the proton.^ 

VIII. CONCLUSION 

In this article, we have described the nonperturbative QCD vacuum in terms of g-theory, 
where q is identified with the particular gluon condensate (12.41) . A crucial role is played by 



the QCD trace anomaly [13], whose potential relevance to the cosmo 



ogical constant problem 



151 ] (see also Ref. [25| for a discussion in the 



has previously been emphasized in, e.g., Ref. 
context of QED). 

The static equilibrium g-theory gives a gravitating vacuum energy density which is ex- 
actly zero, pvac = 0, according to (16. ip . But in a nonstatic situation (e.g., that of the 
expanding Universe), the gluon condensate is perturbed and a nonzero gravitating vacuum 
energy density results, Pvac 7^ 0. The theoretical value for the remnant vacuum energy den- 



^ After the completion of an earlier version of the present article, we became aware of recent work 23 1 
on the evolution of an expanding FRW universe with a hypothetical decaying vacuum energy density 
proportiona]Ll£) the Hubble parameter. 

^ It is known [24] that the lifetime of a main-sequence star is, to a large extent, determined by the proton 
mass and Newton's constant: t^, ^ -E-pianck/^^p' ^^^^ an additional numerical factor of the order of 10'* 
on the right-hand side for a solar-mass star. This numerical factor contains, however, the fine-structure 
constant a and the proton-electron mass ratio mp/me, so that non-QCD physics enters the estimate. The 
estimate for t^, is, therefore, suggestive but inconclusive from our point of view. 



13 



sity (cosmological constant) is estimated to be given by (16.71) and appears to be of the order 



of the experimental value 
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2l|. 



However, a reliable calculation of the present vacuum energy density p^ac will require a 
detailed study of the gluon-condensate dynamics in an expanding universe. Even though 
that study has barely started and many questions remain, it is remarkable and encouraging 
that an explanation of the so-called "dark energy" can perhaps be found in known physics, 
classical general relativity and quantum chromo dynamics. 
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APPENDIX A: REMNANT VACUUM ENERGY DENSITY FROM QCD 

In this appendix, we give a heuristic derivation of expression (16. 3p for the remnant vacuum 
energy density from the gluon condensate of quantum chromodynamics in an expanding uni- 
verse. The main idea is that nonanalytic behavior of the gravitating vacuum energy density 
as a function of the Hubble parameter H{t) = [da/dt)/a may come from the singularity of 
the gluon propagator in the infrared. 

In order to see how this may happen, it is convenient to use the Gribov picture of 



confinement 
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28|. In the Coulomb gauge, the effective gluon mass would then depend 



on the three- momentum k and would increase in the infrared limit /c = |k| — as 

m(k)~A2QCDA. (Al) 

Such a momentum-dependent mass would come from the instantaneous Coulomb interaction 
between the color charges of the gluons, i.e., from the interaction potential f/(r) oc 1/r in 
coordinate space or [/^ oc in momentum space. 

Here, we consider the possible effects from a more singular behavior of m(k) at extremely 
small k. Assuming a linear confinement potential between gluons U{y) oc Aq^^j^t, with 
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Fourier transform [/k oc A.^Qj^/k'^^ we have the following behavior of the effective gluon mass 
in the extreme infrared region: 

m(k) ~ K%^^/e . (A2) 

I — I 

Recall that, on the one hand, the Richardson potential [29] with a behavior of the 
effective gluon mass gives a reasonable description of heavy-quark systems and that, on 



the other hand, a large- A^^c master field has been suggested ll|], which gives this very same 



potential for color sources in arbitrary nontrivial representations of SU{Nc). Current lattice- 



gauge-theory simulations [28| appear to support the behavior flAip . but are by necessity 
limited to rather small volumes of the order of (1 fm)'^. The conjectured behavior (]A2p 
would hold over length scales L larger than 1 fm (perhaps L > 10 fm) and may provide an 
incentive to push the pure-gauge lattice simulations to their limit. 

In the cosmo logical context, a natural infrared cutoff for the divergent gluon mass is 
provided by the Hubble expansion, 

m{k,H)r^Alc^/{P + H'). (A3) 

The contribution of the Hubble expansion to the vacuum energy density can be estimated 
by using, for example, the zero-point energy of the gluon field. For the FRW universe (or, 
more specifically, the de-Sitter universe), the estimated contribution of zero-point energies 
from (1A3P is 

pUH) ~ J ^ (m(k, H) - m(k, 0)) 1— \H\ A^^d , (A4) 

where the factor — 1 counts the number of gluons in a pure SU{Nc) Yang-Mills theory. 
As argued in the main text, the vanishing of the gravitating vacuum energy density Pvac in 
Minkowski spacetime {H = 0) would be due to the self-adjustment of a g-type variable. 



Even though flA4|) has the wrong sign (cosmology [20|, l2l|, |2^ suggests pvac = —-Pvac > 0), 
the important point is to have found that a term of order \H\ Aq^j-, can arise at all. The 
contributions of the fermionic quarks, which have not been considered up till now, may, 
in principle, reverse the overall sign of flA4l) . In any case, the zero-point-energy estimates, 
which are applicable to equilibrium vacua, have only heuristic value if the dynamics of 
the nonequilibrium vacuum is considered: these estimates may give the correct order of 
magnitude but not the exact number or even the sign. 
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Turning to the dynamics, the infrared behavior of QCD in (lA2p induces nonanalytic 
higher-order derivative terms in the gradient expansion mentioned in Sec. IIV[ The singu- 
lar infrared behavior also leads to nonanalytic higher-derivative terms in the master-field 
equation relevant to right-hand side of (14 .Sap . This could, for example, give^ 

G^^D^G^" = ci TQCD^D. (g^, G^'^) + ■ • • , (A5) 

with a numerical coefficient ci, the microscopic time scale tqcd ~ and the invariant 

d'Alembertian □ defined in terms of the master gauge field A^{x) and the standard affine 
connection r^^(x) from the metric g^u{x) and its inverse 

For a fiat FRW universe with a time-dependent homogenous master field fl4.1bl) and 
corresponding scalar field Ji{t) from (]4.8b[) . the differential equation (14 .Sal) can then be 
approximated as 

djlit) 



dt 



-Ajlit) 7(t) H{t)' TQCD + 0{H' t^cd) , (A6) 



with a factor \H\ tqcd in the first term on the right-hand side from the nonanalytic higher- 
derivative term in (]A5|) and a dimensionless function ^{t) from the full master-field dynamics. 
For comparison, analytic higher-derivative terms would give the much smaller factor Tqcd 
contained in the second term on the right-hand side of flA6l) . 

The present Universe at coordinate time t = to > (setting the Big Bang coordinate 
time to zero, teB = 0) has a Hubble constant Hq = //(to) ~ 1/to > and may be considered 
to have a vacuum state near equilibrium, Jt = fiQ + 6fi, for \6fi/fio\ <^ 1 and fiQ given by 
(16. lap . The ordinary differential equation (1A6P gives then approximately 

5ii pa -4/io 7(to) to tqcd ~ -/io 7o \Ho\ tqcd , (A7) 
with all numerical factors absorbed in the constant 70 and using the positivity of to- The 



^ The quartic root of the differential operator V = □ ^ in (|A5ll can be defined as V^^'^ = 
hm^io (2V2/7r) V dk (fc^+P) ^ . Note that the eigenvalues of I? are non- negative, also for a spacetime 
metric with Lorentzian signature. 
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chemical-potential shift (1A7P results in the following nonzero vacuum energy density dOH): 
Pvac(/i) ~ eo - &i g(/io) - hi {dq/d^) 6fi 
~ (^0 Xo) 

~ ^1 (^qcd/^i) (mol^ol/AQCo) 



7o /^o Aqcd |ifo 



^ 7o6i A^cdI^oI, (A8) 

where the derivative of (]4.8bp with respect to q has been used in the second step, the 
combined results (16. Id) . (15.41) . (16.21) . and flA7l) in the third step, and (I6.1al) in the last step. 
The final expression flASp . with positive 7o &i for a de-Sitter-like universe (ifdcs ~ -f^o > 0), 
is precisely of the form (16. 3p with / = 7o &i. Purely theoretically, the first two steps in (lASP 
highlight the importance of the vacuum compressibility xo = xiQo) fo^^ the dynamics of the 
vacuum energy density, which has also been noted, for a different model, in Eq. (2.9) of 
Ref. [l9|. 

Result (lASp or equivalently (16.30 corresponds to a nonanalytic term in phenomeno- 

logical f{R) modified-gravity theories, where a tilde has been added to the function f{R) in 
order to distinguish it from the numerical factor / used elsewhere in this article (see Ref. 
for references on this type of modified-gravity theories). Specifically, the f{R) gravity in- 
duced by QCD is given by 

f{R) = -R-M ^/\R\ + ■ ■ ■ , (A9) 



with M > and the same conventions for the Ricci scalar R as in Refs. 3, 3, 4]. The 
term in flA9p stands for all terms with G Z, while the ellipsis indicates 

other higher-order terms in R. Note that the complete gravitational action may also have 
particular terms involving the Ricci tensor and Riemann tensor, which, for the de-Sitter 
metric, give a vanishing contribution to the generalized Einstein field equation. 

The modified-gravity model with f{R) from (1A9P belongs to the class of chameleon-type 



models 



30 
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32j. For the case of the gluon-condensate vacuum, the corresponding mass 



scale M in (IA9p is given by 

M ~ / A?,co/B?,„. - 2 X 10-»eV (^) ) , (AlO) 
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which, up to a factor An, corresponds to the inverse of (17. ID . At small curvatures, \R\ < M^, 
the square-root term in (\A9\i becomes significant and leads to a large-distance modification 
of gravity due to the existence of a gluon condensate. It will be of interest to work out the 
details of the corresponding cosmological model. 



[1] F.R. Klinkhamer and G.E. Volovik, "Self-tuning vacuum variable and cosmological constant," 

Phys. Rev. D 77, 085015 (2008). [a?Xrv:0711.31701 
[2] F.R. Klinkhamer and G.E. Volovik, "Dynamic vacuum variable and equilibrium approach in 

cosmology" Phys. Rev. D 78, 063528 (2008), larXiv:0806^2805J 
[3] F.R. Khnkhamer and G.E. Volovik, "/(i?) cosmology from g-theory" JETP Lett. 88, 289 

(2008), larXiv:0807.38Ml 
[4] S. Weinberg, Cosmology (Oxford University Press, Oxford, England, 2008). 
[5] T.P. Cheng and L.F. Li, Gauge Theory of Elementary Particle Physics (Clarendon Press, 

Oxford, England, 1984). 

[6] (a) M.A. Shifman, A.I. Vainshtein, and V.I. Zakharov, "QCD and resonance physics: Theo- 
retical Foundations," Nucl. Phys. B 147, 385 (1979); (b) M.A. Shifman, A.I. Vainshtein, and 
V.I. Zakharov, "QCD and resonance physics: Apphcations," Nucl. Phys. B 147, 448 (1979); 
(c) M.A. Shifman, Vacuum Structure and QCD Sum Rules (North-Holland Publishing Co., 
Amsterdam, The Netherlands, 1992). 

[7] S. Elitzur, "Impossibility of spontaneously breaking local symmetries," Phys. Rev. D 12, 3978 
(1975). 

[8] G.K. Savvidy, "Infrared instability of the vacuum state of gauge theories and asymptotic 
freedom," Phys. Lett. B 71, 133 (1977). 

[9] M. Novello, S.E. Perez Bergliaffa, and J. Salim, "Nonlinear electrodynamics and the acceler- 
ation of the Universe," Phys. Rev. D 69, 127301 (2004), arXiv:astro-'ph703 120 93] 
[10] D.N. Vollick, "Homogeneous and isotropic cosmologies with nonlinear electromagnetic radia- 



18 



tion," Phys. Rev. D 78, 063524 (2008). [MXrv:0807.0448[ 
[11] J. Greensite and M.B. Halpern, "Suppression of color screening at large A^," Phys. Rev. D 
27, 2545 (1983). 

[12] D.J. Gross and F. Wilczek, "Ultraviolet behavior of non-Abelian gauge theories," Phys. Rev. 

Lett. 30, 1343 (1973); H.D. Politzer, "Reliable perturbative results for strong interactions?" 

Phys. Rev. Lett. 30, 1346 (1973). 
[13] N.K. Nielsen, "The energy-momentum tensor in a non-Abelian quark gluon theory," Nucl. 

Phys. B 120, 212 (1977); J.C. Collins, A. Duncan, and S.D. Joglekar, "Trace and dilatation 

anomalies in gauge theories," Phys. Rev. D 16, 438 (1977). 
[14] A.R. Zhitnitsky, "Vacuum energy, EoS, and the gluon condensate at finite baryon density in 



QCD," AIP Conf. Proc. 892, 518 (2007), arXiv:hep-ph/0701065 



[15] R. Schiitzhold, "Small cosmological constant from the QCD trace anomaly?" Phys. Rev. Lett. 
89, 081302 (2002). 

[16] J.D. Bjorken, "The classification of universes," |arXiv:astro-ph/0404233. 

[17] G.S. Bah, "QCD forces and heavy quark bound states," Phys. Rept. 343, 1 (2001), 



arXiv:hep-ph/0001312 



[18] Y.B. Zeldovich, "Cosmological constant and elementary particles," JETP Lett. 6, 316 (1967). 
[19] F.R. Klinkhamer, "Equilibrium boundary conditions, dynamic vacuum energy, and the big 



bang," Phys. Rev. D 78, 083533 (2008), arXiv:0803.0281[ 

[20] E. Komatsu et al, "Five-year Wilkinson Microwave Anisotropy Probe (WMAP) observations: 
Cosmological interpretation," Astrophys. J. Suppl. 180, 330 (2009). larXivrOSOgTOSiTl 

[21] A. Vikhlinin et al, "Chandra Cluster Cosmology Project III: Cosmological parameter con- 
straints," Astrophys. J. 692, 1060 (2009) larXiv:0812.2720l 

[22] A.G. Riess et a/., "New Hubble Space Telescope discoveries of type la supernovae at z > 1: 
Narrowing constraints on the early behavior of dark energy," Astrophys. J. 659, 98 (2007), 



arXiv:astro-ph/0611572 



[23] S. Carneiro, C. Pigozzo, H.A. Borges, and J.S. Alcaniz, "Supernova constraints on decaying 



vacuum cosmology," Phys. Rev. D 74, 023532 (2006), arXiv:astro-ph/0605607 ; H.A. Borges, S. 
Carneiro, J.C. Fabris, and C. Pigozzo, "Evolution of density perturbations in decaying vacuum 



19 



cosmology," Phys. Rev. D 77, 043513 (2008), larXiv: 071 1.2 689': H.A. Borges, S. Carneiro, 
and J.C. Fabris, "Evolution of density perturbations in decaying vacuum cosmology: The 
case of nonzero perturbations in the cosmological term," Phys. Rev. D 78, 123522 (2008), 
larXiv:0809.4939[ 

[24] R.H. Dicke, "Dirac's cosmology and Mach's principle," Nature 192, 440 (1961); J.D. Barrow 
and F.J. Tipler, The Anthropic Cosmological Principle (Clarendon Press, Oxford, England, 
1986), Sees. 4.6 and 5.6. 

[25] L. Labun and J. Rafelski "Trace anomaly of nonlinear electrodynamics and its 

(anti)gravitational effect," larXiv:081 1.44671 
[26] V.N. Gribov, "Quantization of non-Abehan gauge theories," Nucl. Phys. B 139, 1 (1978). 
[27] D. Zwanziger, "Lattice Coulomb Hamiltonian and static color-Coulomb field," Nucl. Phys. B 



485, 185 (1997), arXiv:hep-th/9603203 



[28] G. Burgio, M. Quandt, and H. Reinhardt, "Coulomb gauge gluon propagator and the Gribov 

formula," Phys. Rev. Lett. 102, 032002 (2009), [ a?Xrv :0807.32911 
[29] J.L. Richardson, "The heavy quark potential and the T, J/^ systems," Phys. Lett. B 82, 272 

(1979). 

[30] J. Khoury and A. Weltman, "Chameleon cosmology," Phys. Rev. D 69, 044026 (2004), 



arXiv:astro-ph/0309411 



[31] T. Faulkner, M. Tegmark, E.F. Bunn, and Y. Mao, "Constraining f{R) gravity as a scalar- 



tensor theory," Phys. Rev. D 76, 063505 (2007), j arXiv: astro-ph/0612569 [ 
[32] P. Brax, C. van de Bruck, A.C. Davis, and D.J. Shaw, "/(i?) gravity and chameleon theories," 
Phys. Rev. D 78, 104021 f2008^. [a?Xrv:0806.3415[ 



20 



